On 
ON 



X 



MAD/PH/822 

University of Wisconsin - Madison YUMS 94-11 

SNUTP 94-34 

(July 1994) 



^- ■ Testing J/ip Production and Decay Properties 



in Hadronic Collisions 

C.S.Kim a and E. Mirkes 6 
Department of Physics, Yonsei University, Seoul 120-749, KOREA 



£> ' b Physics Department, University of Wisconsin, Madison WI 53706, USA 

00 ' 

m 

o 

■"sj- ■ The polar and azimuthal angular distributions for the lepton pair arising from the 

On 

decay of a J/tp meson produced at transverse momentum p T balanced by a photon 

^ . [or gluon] in hadronic collisions are calculated in the color singlet model (CSM). It 

^S^ . is shown that the general structure of the decay lepton distribution is controlled by 



four invariant structure functions, which are functions of the transverse momentum 
and the rapidity of the J/ip. We found that two of these structure functions [the 
longitudinal and transverse interference structure functions] are identical in the CSM. 
Analytical and numerical results are given in the Collins-Soper and in the Gottfried- 
Jackson frame. We present a Monte Carlo study of the effect of acceptance cuts 
applied to the leptons and the photon for J/ip + 7 production at the Tevatron. 



I. INTRODUCTION 

The study of the properties of the bound states of heavy quarks plays a central role 
in the understanding of perturbative Quantum Chromodynamics. Besides of the intrinsic 
interest in the bound state properties of a J/0, the study of the 3 S\ — cc bound state is 
also of particular importance in the measurement of b quark production. So far there has 
been intensive experimental studies of the J/0 production rate and transverse momentum 
distributions in hadronic collisions both at UAl P] and CDF 0. However, the observed 
J/0 rate was found to be markedly higher than the predicted one's in [[3], where in addition 
to the direct charmonium production also contributions from the x an d the production 
resulting from B decays were taken into account. It has recently been pointed out in 0] 
that at large transverse momentum of the J/0 an additional production mechanism comes 
from fragmentation contributions of a gluon or a charm quark into charmonium states. 

In this paper we propose to measure the decay lepton distribution of the J/0 as a detailed 
test of the production mechanism of the cc bound state. We restrict ourself in the present 
paper to the study of J/ip produced in association with a photon. This has several advantages. 
First of all the experimental signature of a J/0 [decaying into an e~e + or \f~ fi + pair] and 
a 7 with balancing transverse momentum is a very clean final state. Second, there are no 
fragmentation contributions or contributions from radiative x decaysa to J/ip + 7 production 
and the only subprocess contributing to a J/0+7 final state is the gluon gluon fusion process. 

For J/0's produced with transverse momentum p T [balanced by the additional photon] 
one can define an event plane spanned by the beam and the J/0 momentum direction which 
provides a reference plane for a detailed study of angular correlations. The decay lepton 
distribution in J/0 — > l~l + in the J/0 rest frame is determined by the polarization of the 
J/0. Therefore, the study of the angular distribution can be used as an analyzer of the 



x The radiative \J decays can produce J/0 at both low and high p T , but the photon produced 
will be soft [E ~ 0(400 MeV)]. 



J/0-polarization. It is thus possible to test the underlying J/^-pro duct ion dynamics [in 
our case the color singlet model (CSM) ||] in much more detail than is possible by rate 
measurements alone. It is shown that the angular distribution of the leptons in the J/0-rest 
frame has the general form 

a ~ (1 + cos 9 2 ) + - A (1 - 3 cos 2 9) + A 1 sin 29 cos + - A 2 sin 2 9 cos 20 , (1) 

where 9 and denote the polar and azimuthal angle of the decay leptons in the J/ip-rest 
frame. The coefficients Ai in Eq. ([!]) are functions of the transverse momentum and rapidity 
of the J/ip. We present analytical results for these coefficients in the Collins Soper (CS) 
H and Gottfried Jackson (GJ) frame. The functions A and A 2 are found to be equal in 
the CSM. An experimental analyses of this equality would be a sensitive test of the CSM. 
We show numerical results for the transverse momentum [pt(J/0)] an d rapidity [y(J/ip)} 
dependence of the coefficients Ai and study the effect of the acceptance cuts on the and 
cos# distributions in the CS and GJ frame. 

The coefficients Ai are ratios of helicity cross sections [see Eq. (f25|)] and the uncertainties 
from higher order QCD corrections and the wave function of the bound state are expected 
to cancel in this ratio. This has been shown for the related process of high px lepton 
pair production in the Drell-Yan process [0] , where large QCD corrections to the individual 
helicity cross sections almost cancel in the corresponding ratios Af Y . A similar effect is 
expected in the case of J/ifj production and our LO predictions for the Ai and the [normalized] 
and cos 9 distributions provide therefore a useful test of J/0 production. 

Our analytical results for the coefficients Ai are also valid for J/ip + g production within 
the CSM. As mentioned before, there are additional contributions to this final state |4[]. A 
study of the angular distribution as proposed in this paper may allow to disentangle the 
different contributions. 

Recently, J ftp's produced in association with a photon has been proposed as a clean probe 
to study the gluon distribution and the polarized gluon distribution of proton as well as to 
investigate qualitatively the production mechanism of heavy quark bound-states [f|. An 



investigation of the decay lepton distribution and structure functions of J fa's production 
and decay in e~e + -experiments is presented in ||. 

The remainder of this paper is organized as follows: In Sec. II, analytical results for the 
helicity cross sections and the coefficients A{ are derived in the CSM for the CS and GJ 
frame and the formalism for describing the angular distributions is discussed. Numerical 
results are presented in Sec. Ill and summary remarks are given in Sec. IV. Some technical 
details of the calculation are relegated to the appendix. 

II. STRUCTURE FUNCTIONS AND ANGULAR DISTRIBUTIONS 

We consider the angular distribution of the leptons coming from the leptonic decay of 
J/0's produced with non-zero transverse momentum in association with a photon in high 
energy proton-antiproton collisions. For defmiteness we take 

p(Pi) + p(P 2 ) - Jft{P) + l(k) +X^r(l) + /+(/') + 7 (A;) + X , (2) 

where the quantities in the parentheses denote the four-moment a. In leading order of pertur- 
bative QCD [O(a^)], J/ip + 7 can only be produced in gg fusion: 

9iPx) + 9(P>) - JH>{P) + l{k) +X. (3) 

In the parton model the hadronic cross section is obtained by folding the hard parton 
level cross section with the respective parton densities: 

dfidydn* = I ^W^i,/4)/ 2 (*2,/4) ^^ (*iPi,*Aa.(/4)) , (4) 

where g h (x,fi 2 F ) is the probability density to find a gluon with fraction x in hadron h if it 
is probed at a scale p? F . The partonic cross section for the process in Eq. (||) is denoted by 
da. Denoting hadron level and parton level quantities by upper and lower case characters, 
respectively, the hadron and parton level Mandelstam variables are defined by 

S=(P 1 + P 2 ) 2 , T=(P 1 -P) 2 , U=(P 2 -P) 2 , (5) 



and 

s = (p 1 +p 2 ) 2 = x 1 x 2 S , 

t = ( Pl -P) 2 = Xl (T-P 2 ) + P 2 , (6) 

u = (P 2 --P) 2 = x 2 (U-P 2 ) + P 2 , 

where p\ = X\P\ and p 2 = X\P 2 . The rapidity y of the J/ip in the laboratory frame can be 
written as 

1 , ( p2 ~ U \ 
y= 2 g \ P 2 -t ) ' ^ 

and the transverse momentum of the J/ip is related to the Mandelstam variables via 

(P 2 -U)(P 2 -T) p2 _ut 
Pt S s ' K ' 



The angles # and <fi in dO* = dcos8d(fi in Eq. (|) are the polar and azimuthal decay 
angles of the leptons in the J/ip rest frame with respect to a coordinate system described 
below. The angular dependence in (||) can be extracted by introducing four helicity cross 
sections corresponding to the four none zero combinations of polarization density matrix 
elements [see Appendix A of this paper, and Appendix C of for details]: 

167r da da u+L . 2 . da L . 2 . . . 

3 dp 2 dy d cos 9 d<p dp 2 T dy dp 2 ^ dy 

+ , - , 2 sin 2 9 cos 2(p + -—— 2\/2 sin 29 cos <p . 
dp z T dy dp z T dy 

The hadronic helicity cross sections d d % a d are again obtained by convoluting the partonic 
helicity cross sections with the gluon densities. One has 

J dx 1 dx 2 g h '{x uf i 2 F )g h ^x 2 ,fi 2 F )'^-. (10) 



dp\dy 

Each of the partonic helicity cross sections is calculated in the CSM. The unpolarized dif- 
ferential production cross section is denoted by a u+L whereas a L)T ' 1 characterize the polar- 
ization of the J/ip, i.e. the cross section for the longitudinal polarized J/0's is denoted by 



a L , the transverse- longitudinal interference by a 1 , and the transverse interference by a T (all 
with respect to the z-axis of the chosen Jfip rest frame). The results for the helicity cross 
sections a a are dependent of the choice of the z axis in the rest frame of the J/0. We will 
derive explicit results for the Collins-Soper (CS) and the Gottfried- Jackson (GJ) frame. The 
Lorentz transformation matrix for the boost from the laboratory frame to the CS and GJ 
frame is given in Appendix B. 

In the CS frame || the z-axis bisects the angle between Pi and — P 2 

CS : Pi = Pi (sin7cs,0,cos7cs) , 

P 2 = P 2 (sin 7 CS , 0, - cos j C s) , 



(11) 



with 



COS 7c-5 



™%S 



\(T-ml)(U-ml)--\ml+pl 



m l 



(12) 



sin7cs = -yl - cos 2 7c5 • (13) 

In the GJ frame (also known as t-channel helicity frame) the z-axis is chosen parallel to 
the beam axis 



with 



GJ : Pi = Pi (0, 0, 1) , 

P 2 = P 2 (sin7 GJ , 0, cos7 G j) 



2mlS pi ™ 2 



(14) 



C0S7GJ = ' ' (r-,4)(c/-,4) - ^4 • (15) 



sin7 GJ = -yl -cos 2 7 GJ . (16) 

The beam and target energies in the rest frame are 

Pi = (mj - P)/(2m v> ), P 2 = (ml - U)/{2m il ) . 

The partonic helicity cross sections in Eq. (|IUD are calculated applying the technique 
described in [0 to the CSM [see Appendix A for details]. They are given by: 



s °cs,gj 167raafofy 2 2 

^*T- 27, imi H cs,ajS(s + t + u-m,,), (17) 

with 

rrt/+L _ f , * | " /ion 

cs ~ (t-m^) 2 (w-m^) 2 (u-m$) 2 (a-m$) 2 (s-m$) 2 (*-rr4) 2 ' ^ ^ 

l ut ( 4m 2 s 3 c/+l\ 

^ cs = 2(t-mJ)(«-mJ) ^( S -mJ)2(t-mJ)2( u -mJ)2 + ^ cs J ' (19) 

H T CS = \nt , (20) 

j m^y/sut s(s 2 - ut)(t - u) 

V2(s - m\) 2 (t - mj) 3 (tt - m^) 3 ' 

and for the GJ frame 
H U G j L = H%$ L , (22) 



2m 2 b stu(s 2 + -u 2 ) 



^Gi = ^GJ , (24) 

7 m v ,\Zgtw(g-M)[g 2 (M-t) + u 2 (s-t)} 

GJ ~ y/2{s-m%)\t-m%Yiu-m%f ' l j 

As mentioned before, H u+L denotes the matrix element contribution for the production rate 



and was first calculated in [III]. All other matrix elements correspond to the production of 
polarized J ftp's and are given here for the first time. Replacing 16a by 15a s in Eq. (|I7|), the 
results in Eqs. ( |T7| - |2l)D are also valid for J/ip + g production within the CSM. 

Note, that the matrix elements can be expressed in terms of p 2 , y by noting that 

t = m\- \j{m\+p 2 T )S x x e 



-y 

> 



u = m\- J (ml + p 2 )S x 2 e y , (26) 



so that e.g. 



(fn T\\Tn ; ?y1 lit 

— = P 2 T + rri^ , and also — = p 2 T ■ (27) 



From the explicit expressions in Eqs. (|18"l - f2"5D one observes, that the longitudinal and the 
transverse interference cross sections a L and a T are related by <Jl = 2a? and that they vanish 
for p T — > with p 2 . The longitudinal-transverse interference cross section a 1 vanishes with 
p T , whereas the production cross section a u+L for unpolarized J/0's is finite in the p T — > 
limit. 

Introducing standard angular coefficients |J 



An 



2 da L 
da u + L 



A, 



2^ da 1 
da u + L 



A, 



Ada 1 



da u+L 



(28) 



the angular distribution of Eq. (Rf) is conveniently written as: 



da 



3 da u+L 



dp 2 dyd cos 9 d<f> 16ir dp 2 dy 



(l + cos 2 #)+ -A (l-3cos 2 #) 



1 



+ A\ sin26 l cos0 + -A 2 sin 2 9 cos 20 



Integrating the angular distribution in Eq. ( p9|) over the azimuthal angle yields: 



(29) 



dp?, dy d cos 9 



C (1 + a cos 



(30) 



where 



C 



3 da c/+L 



An 

1 + — 
2 



8 dp 2 , dy 
And integration over ^ yields 

da 1 da c/+L 



a: 



2-3A 

2 + An 



(31) 



where 



dp 2 , dy d(f) 2tt dp 2 dy 



?= A T- 



l+/3cos20) 



(32) 



(33) 



Before discussing numerical results, let us briefly discuss a possibility strategy for extracting 
angular coefficients. By taking moments with respect to an appropriate product of trigono- 
metric functions it is possible to disentangle the coefficients A\. A convenient definition of 
the moments through 



8 



_ / do-(pr, y, 9,(p) m d cos d0 
lm; ~ fda(p T ,y,0,<f>) d cos dd(f> ' l j 



which leads to the following results 



(1) = 1, (35) 



.1,.. „ 2,,, 3 r. 2 



(-(l-3cos 2 0)> = - (A- -J , (36) 

(sin 20 cos 0) = - At, (37) 



5 

1 

lb 



(sin 2 cos 20) = — A 2 . (38) 



III. NUMERICAL RESULTS 

We will now present numerical results for J/0 + 7 production at the Tevatron collider 
center of mass energy [\^S = 1.8 TeV] including the decay J/0 — > ^~yu + . All results are 
obtained using the gluon density parametrization from GRV ]TT| with Ai^L = 200 MeV and 
the one-loop formula for a s with 4 active flavours. If not stated otherwise, we identify the 
renormalization scale fi R and the factorization scale /j, f in Eq. (^) and set them equal to 
ji 2 = H F = H R = (m 2 +p 2 r (J/4>))- The value for the bound state wave function at the origin 
|i?(0)| 2 is determined from the leptonic decay width of J/0: T(J/0 — >■ e~e + ) = 4.72 keV, 
therefore |i?(0)| 2 = 0.48 GeV 3 . 

Fig. 1 shows the transverse momentum distribution of the the J/0 where we have ap- 
plied the cuts I y (7, l~, l + )\ < 2.5 and pt(1~ ,l + ) > 1.8 GeV on the final state particles. The 
p T distributions are leading order predictions and are therefore strongly dependent on the 
choice of the scales /x 2 and the radial wave function of the bound state. In fact, the NMC 



collaboration found |12] that the LO predictions describes the shape of all kinematical vari- 
ables quite well, while the predictions for the normalization of the signal was too small by a 
factor 2-3. This large "K- factor" is probably due to the nonrelativistic treatment of the J/0 
and due to higher order QCD corrections in the CSM. One can expect a similar K-factor for 



9 



our reaction, so that our result of Fig. 1 should be considered as conservative estimates. To 
give a feeling for the scale dependence we show results for [j? — \i 2 F — n 2 R — (m 2 + p 2 T (J/ip)) 
[solid line], fi 2 = (mi + p 2 (J/fy>))/4 [dashed line], and /i 2 = 4(m 2 + p 2 ,(J / V>)) [dotted line]. 

Fig. 2 shows the rapidity distribution of the J/ip normalized to the total production 
cross section. The following cuts are applied: pt(j,1~ ,l + ) > 3 GeV. The shape of this 
distribution is quite sensitive to the gluon density function, and therefore useful to extract 
the gluon density at the small x f|. 

In Fig. 3, we show the transverse momentum distributions of the decay leptons from the 
J/ip after the cuts: \y(jJ-,l + )\ < 2.5, p T (l~,l + ) > 1.8 GeV and p T {l) > 2 GeV. We also 
plot the spectrum of the harder (denoted by 'b') and softer (denoted by V) lepton separately, 
where we did not distinguish between the charge of the leptons. Note that together with 
the hard photon, this harder lepton (denoted by 'b') can be used to construct a trigger for 
our reaction. Like in Fig. 1, the pt distributions of the leptons are strongly dependent on 
the choice of the scales fi 2 and the radial wave function of the bound state. 

In Figs. 4 and 5 we show numerical results for the coefficients Aj in Eq. ( |2"§D as a function 
of Pt(J/iP) and y(J/ip) in the CS [Figs. 4, 5(a)] and GJ [Figs. 4, 5(b)] frame. These coefficients 
have been extracted from the Monte Carlo program by using the moments defined in Eq. (|J3|) . 
No acceptance cuts have been applied to the leptons or the photon. One observes that the 
coefficients are strongly dependent on priJ/^) both in the CS and GJ frame. As mentioned 
before, the coefficients A and A 2 are exactly equal in lowest order in both lepton pair rest 
frames. The angular coefficient Ai is zero in the CS frame for all values of p T (J/ip). The 
reason is that the matrix element for A% is antisymmetric in u and t [see Eqs. (^TJ,^8|)] and 
therefore in x\ and x 2 , whereas the product of the gluon distributions is symmetric under 
the interchange of x\ and x 2 . As a consequence, the rapidity distribution for Ai in the CS 
frame [Fig. 5(a)] is also antisymmetric around y = 0. However, this is different for the GJ 
frame [see Eq. (p5])]. All coefficients Ai vanish in the limit p T (J/ip) —>■ 0, which can be 
directly seen from our analytical expressions in Eqs. ( P^ - |25]J2S| ) . 

A similar relation A® Y = A 2 Y was found in LO [0(a s )] in the Drell-Yan process |13| 
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p + p — > V + X — > l + l' + X, where V denotes a gauge boson produced at high pt- In [[?], 
the complete NLO corrections to the coefficients Af Y are calculated and the corrections 
are found to be fairly small for the ratio's Ai of the corresponding helicity cross section. 
We expect also here, that the LO results for the ratio's Ai in J/0 production are almost 
not affected by higher order QCD corrections. Note also, that the coefficients Ai are not 
dependent on the bound state wave function. The measurement of these coefficients would 
be a sensitive test of the production mechanism of J/0's. However, as we will see later, 
experimental cuts introduce additional complicated angular effects and the resulting data 
sample can no longer be described by the simple angular distribution in Eq. (f29|) . 

Figures 6-11 show the normalized and cos# distributions for the leptons in the CS and 
GJ frame for three different bins in the transverse momentum of the J/0. To demonstrate 
the effects of acceptance cuts, results are shown first without cuts and then with typical 
acceptance cuts imposed to the leptons and the photon. 

Figures. 6 and 7 show the and cos 9 distributions in the CS and GJ frame without cuts. 
The curves in Figs. 6 and 7 can be obtained from the results in Fig. 4 using the coefficients 
a and (3 defined in Eqs. ([J(],[3^). One observes a fairly strong dependence both in the CS 
[Fig. 6(a)] and the GJ [Fig. 7(a)] frame for the three p t (J/ip) bins, which is expected from 
(3 = A2/4 and the results in Fig. 4. The dependence is larger for high py(J/0) bins in 
the CS frame, whereas an opposite behavior is found for the GJ frame. The coefficient a in 
Eq. (|30|) is negative for the three bins in the CS frame and positive in the GJ frame. The 
corresponding cos 8 distributions are therefore decreasing [increasing] with cos# in the CS 
(Fig. 6(b)) [GJ (Fig. 7(b))] frame. Let us recall, that the angular distributions in Figs. 6 
and 7 are determined by the polarization of the J/0 with respect to the z-axis in the CS and 
GJ frame [see Eq. (Q]. 

Figs. 8 and 9 show the and cos 9 distributions in the CS and the GJ frame for the same 
bins in py(J/0) as in Figs. 6 and 7 but with the cuts 

p T (r i l + )> 1.8 GeV |y(7,r,Z+)| <2.5. (39) 
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The cuts have a dramatic effect on the shapes of the distribution. The shapes of the dis- 
tributions are now governed by the kinematics of the surviving events. □ The cuts, which 
are applied in the laboratory frame, introduce a strong and cos 9 dependence. The cuts 
remove mainly events around = 0, n, 2i\ [Figs. 8, 9 (a)] and cos 9 > 0.5 for the low pr(<W0 
bin [solid line in Figs. 8, 9(b)]. The distributions in Figs. 8, 9(a) are very different from 
the "dynamical" distribution in Figs. 6, 7(a) and there is only very little sensitivity to the 
polarization dependence (shown in Figs. 6, 7) left in Figs. 8 and 9. 

We have also analyzed the effect of the cuts separately by using the correct matrix element 
for J/ip production, but with isotropic decay of the J/ifi, i.e., neglecting spin correlations 
between J/ip production and decay. The angular distributions in this case are very similar 
to the ones shown in Figs. 8, 9 for the full matrix element. 

In Figs. 10 and 11 we show ratios of the and cos 9 distributions in the CS and GJ 
frame for the same bins in j>r( J/0) as in Figs. 6-9; the distribution with full polarization has 
been divided by the distribution obtained with isotropic decay of the J/ip. Cuts are included 
in both cases. The large effects from the cuts are expected to almost cancel in this ratios. 
In fact, we nearly recover the and cos 9 dependence of Figs. 6 and 7 in Figs. 10 and 11. 
Especially, the two high Pt(J/V0 bins [dashed and dotted lines in Figs. 10 and 11] contain 
most of the polarization dependence seen in Figs. 6 and 7. Comparing the results in the 
CS [Figs. 10 and 6] and the GJ frame [Figs. 11 and 7], the ratios in the CS frame are more 
sensitive to the distribution, whereas the ratios in the GJ frame are more sensitive to the 
cos 9 distribution, when cuts are applied. 

Therefore, to regain sensitivity to the polarization effects in the presence of large kine- 
matic cuts, we propose to divide the experimental distributions by the Monte Carlo distri- 



We thank V. Barger and R. J. N. Phillips for pointing out this fact. 
3 Analogous ratios have been recently applied in [14| to regain the sensitivity to polarization effects 



in W and Z production in hadronic collisions. 
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butions obtained using isotropic J/0 decay. 

IV. SUMMARY 

Hadronic J/0+7 production has been evaluated in the nonrelativistic bound state model. 
In leading order this final state can only be produced by gluon fusion. Analytical formula for 
the decay lepton distributions in terms of four structure functions are presented. If no cuts 
are applied, the angular distribution in the J/0 rest frame is determined by the polarization 
of the J/0. We present Monte Carlo studies of the leptonic decay products of high p? J/0's 
produced in association with a photon, when cuts are imposed on the photon and leptons. 
When acceptance cuts are imposed on the leptons, the angular distributions are dominated 
by kinematical effects rather than polarization effects. Polarization effects can be maximized 
by minimizing the cuts. Alternatively, it may be possible to retain J/0 polarization effects 
by "dividing out" the kinematic effects, i.e., if the histogrammed data are divided by the 
theoretical result for isotropic J/0 decay. 
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APPENDIX A: STRUCTURE FUNCTIONS AND HELICITY CROSS SECTIONS 

In this appendix, we will derive the partonic helicity cross section in Eqs. ( |TT| - p5D for the 
production of polarized J/0's 

The production cross section for 



g(pi)+g(p2)->J/HP)+i(k) 



(Al) 



is 



da = ^-\M\ 2 dPS {2) 



with 



dPS^ 



8tts 



6(s + t + u- P 2 ) dtdu . (A2) 



The amplitude which describes the coupling of the jgg to the J/ip can be calculated within 
the bound state formalism [|[ to be 



j^ggi = _. 



R{0)Tr\G°(^-M t )(-^) 



(A3) 



2 V^ttM^, 

where -R(O) denotes the radial wave function of the bound state, which can be calculated 
either from potential models or can be related to the leptonic decay rate 



|i?(0)p 



M l r 
4a 2 (2/3) 2 ee 



(A4) 



In Eq. ( |A3|) E denotes the J/ip polarization. The amplitude O can be obtained from the 
amplitude for the gluon gluon fusion into a free quark pair at threshold and a photon: 



Q c 



+ 



t^If - 2jji + m*)f{-I? + 2j> 2 + m^ 2 
Ppx Pp 2 
f{-P + 2& + 2j) 2 + m^j-I} - 2j> 2 + m^ 2 

(Ppx + Pp 2 - 2p x p 2 ) Pp 2 
t x {tf - 2^ + m^ 2 {P - 2j> x - 2p\ 2 + m^*' 



Ppi (Ppx + Pp 2 - 2pYp 2 ) 



+ [1 <-> 2] . 



(A5) 
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Here Pi,P2,£i,£2 are the momenta and polarization vectors of the two gluons and e is the 
polarization vector of the photon. Coupling constants, spin average factors and color matri- 
ces contribute a factor (4iia s e2/3) 2 (l/4)(2/3) to production cross section. After summing 
over the polarizations ei a , £2/3, £/* we define the density matrix elements of the J/ip as 

H™' = E s (a)T S5 ' E* 5l (a') , (A6) 



where the hadronic tensor T ss is 



rpSS' = 1 Tr 

" 2048m^ 



O a ^{If - m^ 6 } Tr \O a ^(tf- m^ 5 '} * , (A7) 



and 

E s (±) = 4-(0;±l,-i,0), 

(A8) 

E s (0) =(0; 0,0,1), 
are the polarization vectors for the J/ip defined with respect to some coordinate axis in its 
rest frame (see below). 

The angular dependence of the decay leptons from the J/0 can be extracted by introduc- 
ing the following linear combinations of the density matrix elements [We follow the technique 
presented in Appendix C of [[7]]: 

h u+l = #00 + H ++ + H 



h l = ^00 5 

H T = l(H+- + H~+), 

H 1 = \ (H+° + H°+ - H~° - H°~ .) 



(A9) 



■1 
Collecting all remaining factors, one has 



sda = Ifrraa^ , + ^ rU+l + ^ rL 

dtdu 27 s l 

+2^ sin 29 cos H 1 + 2 sin 2 9 cos 20 H T } 6{s + t + u - P 2 ) . (A10) 



In Eq. ( |A10p 9 and denotes the lepton angles in the lepton pair rest frame. What remains 



is the calculation of the H a in the CS and GJ frame. 
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For this purpose we choose a specific representation for Es(<j)E$,(<j') in Eq. ( JAB] ) and 
express them in terms of of gss',PsPs'jPisPjS' (hj — 1)2). Let us define the following 
covariant projections: 



fP = V$ s ,H ss ', {(3e{U + m u L 2 ,L l2 }) 



(All) 



where 



<r>u+L 

1 55' 



V Ll 

1 55' 



-P L 2 

' 55' 



1 55' 



-955' , 

= -k PuPw , 

= J2 P25p25' , 

1 



(A12) 



-^- 2 (p!5p25' + Pl5'P2s) 



Ex 



P 2 -t 



Eo 



P 2 -u 



(A13) 



We have introduced the hatted tensors 



P5P5' „ p l P 

955' = 955' 7^— , Pi,8 = Pi 5 ]52~Pi $ 



p2 ' -•.- — p2— < ' 2v /p 2 - 2VP 2 

The helicity structure functions H a in Eq. ( |A10| ) for a given J/ip rest frame are linear 
combinations of the covariant projections in Eq. (|A11| ). 
For the CS frame one has [F7fl: 



H u+i 

H L 
H T 
H 1 



\ 



t 







1 





1 

4 cos 2 7cs 
(l+cos 2 7cs) 





1 

4 cos 2 7cs 

(l+cos 2 7 CS ) 







\ 



/ 



-1 



4 cos 2 7cs 
(l-3cos 2 7 CS ) 



2 8 sin 2 "fcs cos2 1CS 8 sin 2 'yes cos2 7CS 8 sin 2 7^3 cos 2 "/cs 

^r-^ , a, . - 1 



CS 



V U i?2 



sin 7cs cos 7cs 



4^ 



sm7cs COS7CS 



H U + L 
H L 12 



\ 



(A14) 



where 



COS 7c5 



m|s 



\ (t ~ m%){u - m%) 



sin 7cs 



1 — cos 2 7cs 



(A15) 
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The explicit results for H^ s are given in Eqs. (|18|-pl|). 
The results for the GJ frame can be obtained from: 

/ }JU+L \ 



and 



H L 
H T 
H 1 





1 



1 (l+cos 2 7 GJ ) 

2 2sin 2 7 G j 



G.I 



V2 



-COS7GJ 
sin jgj 



COS7GJ = 1 



2m\s 



,2\ ' 








1 



CQS7GJ 
sin 2 7g j sin 2 7gj 



2^/2 



sm 7g j 



\ 


( H u+L ^ 




E L X 




H L 2 


/ 


H L 12 
{ J 



(<-mJ)(w-mJ) 
The explicit results for Hqj are given in Eqs. (p^-|25|). 



sm 7 G j 



1 - cos 2 7 GJ 



(A16) 



(A17) 



APPENDIX B: LORENTZ TRANSFORMATIONS TO THE CS AND GJ FRAME 



For completeness, we list here the transformation matrices to the CS and the GJ frame. 
In the laboratory frame, the z-direction is defined by the proton momentum and the re- 
direction is defined by the transverse momentum of the J/ip. The Lorentz transformation 
matrix for the boost from the laboratory frame to the CS and GJ frame is given by [P 
denotes the J/ip momentum] 



' E \ 



Px 

Py 



Pz 

V ) 



( 







p 3 

X.J 1 



Pa 

Vp 1 


P T 


P t Pq 


JC'j' 



pi 









1 



'■' ^ ( E^ 



VP 2 
P T PJ 

Jp^Xt 








CS 



Px 

Py 



Pz 

V ) 



(Bl) 



/«() 
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'e^ 



p,. 



Pz 

V / 



G.I 



where 



Pq 



fp^x% 



/p2 



P?, 



VP~5 


VF 1 




^ 


P T (P0-P3) 


1 





p t (P0-P3) 








1 





p 2 T P0+P 2 P3 


PT 


n 


P 2 P +p|P 3 



fp^x% 



l E \ 



Px 

Py 



Pz 

V ) 



lab 



fB2) 



x t = a /p 2 + ^. 



(B3) 



Figure captions 

Fig. 1 The transverse momentum distribution of the J/0 [or 7] in J/0 + 7 production at the 
Tevatron collider center of mass energy [yS = 1.8 TeV] with subsequent J/0 — ► jjl~ jjl + 
decay. The following cuts have been applied: \y(j, l~, l + )\ < 2.5 and px(l~,l + ) > 1.8 
GeV. The renormalization and factorization scales fi 2 are set equal to (p^il) + m %) 
[solid], 1/4(^(7) +m^) [dashed] and 4 (2^,(7) + m%,) [dotted]. 

Fig. 2 The normalized rapidity distribution of the 7 in J/0 + 7 production at the Tevatron 
collider center of mass energy [v^S = 1.8 TeV] with subsequent J/0 — ► /z~/i + decay. 
The following cuts have been applied: p T (7,Z _ ,Z + ) > 3 GeV. 

Fig. 3 The transverse momentum distribution of the leptons from the decay of the J/0 in 
J/0 + 7 production at the Tevatron collider center of mass energy [yS =1.8 TeV]. 
The following cuts have been applied: |y(7,/",/ + )| < 2.5, pr(l~,l + ) > 1.8 GeV and 
Pt(i) > 2 GeV. Results are shown for the positive (negative) charged lepton [dashed 
lines] as well as the distribution for the lepton with bigger (subscript b) and smaller 
(subscript s) p T without identifying the charge. 

Fig. 4 Angular coefficients Ao, A\ and A^ for J/0 + 7 production and decay in the CS frame 
(a) and GJ frame (b) as a function of the J/0 transverse momentum at \[S =1.8 TeV. 
No cuts have been applied. 

Fig. 5 Same as Fig. 4 for the y(J/ip) distribution. 

Fig. 6 a) Normalized and b) normalized cos 6 distributions of the leptons from J/0 decay 
in the CS frame. Results are shown for three bins in p T (J/0): 
2 GeV < p T (J/*P) < 5 GeV [solid], 
5 GeV < p T (J/ip) < 8 GeV [dashed], 
8 GeV < p T (J/0) [dotted]. 
No cuts have been applied. 
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Fig. 7 Same as Fig. 6 in the GJ frame. 

Fig. 8 Same as Fig. 6 but with the cuts \y(-f,l~ ,l + )\ < 2.5, p T (l~,l + ) > 1.8 GeV. 

Fig. 9 Same as Fig. 7 but with the cuts \y(-f, l~, l + )\ < 2.5, p T (l~,l + ) > 1.8 GeV 

Fig. 10 Ratios of distributions obtained with full polarization effects to those obtained with 
isotropic decay of the J/0. Parts a) and b) are the ratios for the <fi and cos 8 distribu- 
tions in the CS frame, respectively. The cuts \y(j, l~ , l + )\ < 2.5, pt(1~,1 + ) > 1.8 GeV 
are included. 

Fig. 11 Same as Fig. 10 for the GJ frame. 
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